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y—i ' Inclusive bottomonium hadroproduction at the Tevatron is firstly examined in a Monte 

t-H , Carlo framework with the colour-octet mechanism implemented in the event generation. We 

extract some NRQCD colour-octet matrix elements relevant for T(IS') hadroproduction. 
Remarkably we find a quite small contribution (compatible with zero) from feeddown of Xbj 
states produced through the colour-octet mechanism: T(1S) indirect production via XbJ 
decays should be mainly ascribed to the colour-singlet model. Finally we extrapolate to 
O j' LHC energies to predict prompt T(15) production rates. 
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1 Introduction 



B-physics has been reserved an important role in the exciting programme to be developed 
at the Large Hadron Collider or LHC along the 21st century. As an example, the ATLAS 
TDR |]J collects a large number of topics related to charm and beauty flavours allowing 
precise tests of the Standard Model benefitting from a foreseen huge statistics even with 
the machine running at "low" luminosity (~ 10 33 cm -2 

In a series of previous papers @, [5[ |], [5] we examined charmonium hadroproduction 
in a Monte Carlo framework, using PYTHIA 5.7 |J event generator with the colour-octet 
model (COM) J7| implemented in. Basically, such a production mechanism is based on 
the formation of an intermediate coloured state during the hard partonic interaction, 
evolving non-perturbatively into physical heavy resonances in the final state with certain 
probabilities governed by NRQCD This mechanism constitutes a (relativistic) gener- 
alization of the so-called colour-singlet model (CSM) |J which requires the formation of a 
colour-singlet state in the hard interaction itself. Although the discrepancies between the 
CSM and experimental cross sections on bottomonia hadroproduction are smaller than 



those found for charmonia JTO] , still some extra contribution should be invoked to account 
for the surplus observed at the Fermilab Tevatron. 

In this paper we extend our analysis on J/ip and if)' hadroproduction [|| to the bot- 
tomonium family lowest vector resonance, i.e. the T(IS') state. Once again, those matrix 
elements (MEs) determined from Tevatron data in other analysis |Tl|] have to be lowered 
once initial-state radiation of gluons is taken into account. This is because of the raise of 
the {effective) intrinsic momentum (kr) of the interacting partons enhancing the high-p^ 
tail of the differential cross section for heavy quarkonia production (for more details the 
reader is referred to Ref. Q). 

The study of bottomonia production at hadron colliders should permit a stringent 
test of the colour-octet production mechanism, especially regarding the expected (mainly 
transverse) polarization of the resonance created through this process at high-pr- More- 
over, LHC experiments will cover a wider range of transverse momentum than at the 
Tevatron. Therefore, it is worth to estimate, as a first step, the foreseen production rate 
of bottomonium resonances at the LHC and this constitutes one of the goals of this work. 
Thereby any experimental strategy to be conducted in the future (for example a specific 
high-level trigger within the dedicated B-physics data-taking) can be foreseen in advance. 

We have based our study on recent results from Run lb at the Tevatron This 
means significantly more statistics than the data sample from Run la, employed in a 
former analysis [|11|]. However, the different sources of prompt T(IS') production were not 
yet separated along the full accessible pr-range, in contrast to charmonium production. 
Hence we give in this work numerical values for some relevant combinations of long- 
distance MEs (including direct and indirect T(IS') production []) extracted from the fit 
to the CDF experimental points. Nevertheless, we still are able to estimate some colour- 
octet MEs for direct production from the measurements on different production sources 
at p T > 8 GeV 0. 



1 Prompt resonance production includes both direct and indirect channels, the latter exclusively referred to 
feeddown from T(nS) and \bj(nP ) states - i.e. excluding long-lived particle decays 
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2 Implementation of the COM in PYTHIA 



Originally the event generator PYTHIA 5.7 produces direct J /if) and higher Xcj resonances 
via the CSM only ||. It is not difficult to extend this generation to the bottomonium 
family by redefining the resonance mass and wave function parameter accordingly. In our 
analysis we have besides implemented a code in the event generator to account for the 
colour-octet production mechanism via the following a 3 partonic processes for heavy 
quarkonium production: 

9 + g^(QQ)[ 2S+1 Xj] + g (1) 

g + q - (QQ)[ 2S+1 Xj] + q (2) 

q + q^{QQ)\ 2S+1 Xj] + g (3) 

where (QQ) [ 25+1 Xj] stands for a certain heavy quarkonium state denoted by its spectro- 
scopic notation. In particular we have considered the 3 S[ 8 \ x Sq and 3 -Pj 8 ^ contributions as 
leading-order intermediate coloured states. In addition we generated T(IS') and XbjinP) 
(n = 1,2) resonances decaying into Y(IS') according to the CSM as mentioned above. 

A lower p T cut-off was set equal to 1 GeV (by default in PYTHIA) throughout the 
generation since some of the contributing channels are singular at vanishing transverse 
momentum [p. 



2.1 Set of "fixed" and free parameters used in the generation 

Below we list the main parameters, including masses and branching fractions, used in our 
generation with PYTHIA 5.7. We employed the CTEQ2L parton distribution function 
(PDF) in all our analysis. 

Masses and branching fractions: 

• m b = 4.88 GeV 

• rn resonance 2,171^ 

• BR[r(lS)^fi + fi-] = 2.48 % (|I§) 



Colour-singlet parameters (from ||16|| ): 

• < Oj {1S \ 3 Sx) > \ tot = 11.1 GeV 3 , defined as 



< 0? ilS \ 3 S 1 ) > L = £ < 0^ nS) ( 3 S 1 ) > Br[T(nS)^T(lS)X] 

71=1 



• <0^ 1(1P) ( 3 P 1 ) >=6.09 GeV 5 

2 We find from our simulation that gluon-gluon scattering actually stands for the dominant process as expected, 
gluon-quark scattering contributes appreciably however (~ 28% of the colour-octet production cross section), 
whereas quark-antiquark scattering represents ~ 4%. These fractions are slightly larger than those found in 
charmonium production as could be expected from a heavier quark mass 
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• <Of 1{2P \*P x ) >=7.10 GeV 5 

The radial wave functions at the origin (and their derivatives) used in the generation can 
be related to the above matrix elements as 

<0^ 1S \ 3 S 1 )>= ^\R(0)\ 2 (4) 

< of J(nF) ( 3 Pj) > = A( 2 J+ l)|p'(0)| 2 (5) 
whose numerical values were obtained from a Buchmuller-Tye potential model tabulated 



in Ref. [17 



Colour-octet long-distance parameters to be extracted from the fit: 

• < Oj {1S \ 3 Sx) > \ tot , defined as 

< 8 T(ls) ( 3 Si) > \ tot = £{< 0^( 3 S 1 ) > Br[T(nS)^T(lS)X] 

n=l 

+ j:<Or j(nP) ( 3 S 1 )>B r [ Xb j(nP)^T(lS)X]} (6) 
j=o ' 

. < o^es o) > 

. < O^ s \*P ) > 

On the other hand, the differences in shape between the and 3 -Pj 8 ' ) contribu- 

tions were not sufficiently great to justify independent generations for them. In fact, 
temporarily setting < Og 1S \ 3 P ) >= ml < oJ^^^Sq) > and defining the ratio 



da [3p(8)l 
2^ J=Q dpT [ tj 

da [1 5"( 8 )l 



r M = -l (7) 

it is found r ~ 5 as a mean value over the [0, 20] GeV p^-range. Actually the above ratio is 
not steady as a function of the T(15) transverse momentum. Therefore in the generation 
we splitted the pt region into two domains: for px < 6 GeV we set r = 6 whereas for 
Pt > 6 GeV we set r = 4. 

In summary, only the Sq channel was generated but rescaled by the factor r to 
incorporate the 3 Pj 8 ^ contribution as we did in [ffl for charmonium hadroproduction. 



Consequently, in analogy to [11 1 we shall derive a numerical estimate for the combination 



of the colour-octet matrix elements: 



m 2 h 
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2.2 Altarelli-Parisi evolution 



According to the colour-octet model, gluon fragmentation becomes the dominant source 
of heavy quarkonium direct production at high transverse momentum. On the other hand, 
Altarelli-Parisi (AP) evolution of the splitting gluon into (QQ) produces a depletion of 
its momentum and has to be properly taken into account. If not so, the resulting long- 
distance parameter for the 3 S[ 8 ^ channel would be underestimated from the fit [[14] 



The key idea is that AP evolution of the fragmenting gluon is performed from the 
evolution of the gluonic partner of quarkonium in the final state of the production channel 

9 + 9^9*(^(QQ)[ 3 S? ] ]) + 9 (8) 

Let us remark that, in fact, g* is not generated in our code |4j]. Final hadronization 
into a (QQ) bound state is taken into account by means of the colour-octet matrix ele- 



ments multiplying the respective short- distance cross sections [pi] , ||. Nevertheless, it is 
reasonable to assume that, on the average, the virtual g* should evolve at high pt simi- 
larly to the other final-state gluon - which actually is evolved by the PYTHIA machinery. 
We used this fact to simulate the (expected) evolution of the (ungenerated) g* whose 
momentum was assumed to coincide with that of the resonance (neglecting the effect of 
emission/absorption of soft gluons by the intermediate coloured state bleeding off colour 

§)• 

Therefore, event by event we get a correcting factor to be applied to the transverse 
mass of the (QQ) state (for the 3 Sf^ channel only): 



y P T + m \ Q Q) 
(QQ) 

where Pt (Pt) denotes the transverse momentum of the final-state gluon without (with) 
AP evolution and ^(qq) denotes the mass of the resonance. At high pt, 

Pt P = x p x p T (10) 

where pt is the transverse momentum of the resonance as generated by PYTHIA (i.e. 
without AP evolution), whereas for p? < tti^qq) the effect becomes much less significant 
as it should be. Thus the interpolation between low and high pt is smooth with the right 
asymptotic limits at both regimes. 

The above way to implement AP evolution may appear somewhat simple but it remains 
in the spirit of our whole analysis, i.e. using PYTHIA machinery whenever possible. In 
fact, it provides an energy depletion of the fragmenting gluon in accordance with Cho and 
Leibovich's work for charmonium hadropro duct ion jri , [H] . It is worth to note, moreover, 



that the effect of the AP evolution on the generation over the [0,20] GeV pr-range, though 
sizeable, is considerably less pronounced for bottomonium than for charmonium because 
of the larger mass of the former. 

Notice finally that, although we can switch on/off AP evolution and initial-state radi- 
ation at will in the generation, both next-to-leading order effects have to be incorporated 
for a realistic description of the hadronic dynamics of the process. 
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Figure 1: Theoretical curves obtained from a fit using PYTHIA including the colour-octet mech- 
anism for prompt T(1S) production against CDF data at the Tevatron a) without AP evolution 
of the fragmenting gluon, b) with AP evolution of the fragmenting gluon. The CTEQ2L parton 
distribution function and = 4.88 GeV were employed in the fits; dotted line: CSM, dashed 
line: 1 + 3 Pj 8 ^ contribution, dot-dashed line: 3 s[ 8 ^ contribution, solid line: all contributions. 



3 Fit to Tevatron data 



As already mentioned, the theoretical differential cross section on inclusive production 
of prompt T(15')'s stands above Tevatron experimental points for relatively high px if 
the set of long-distance parameters from |]lTJ are blindly employed in the PYTHIA gen- 
eration with initial-state radiation on. Therefore we performed a new fit to recent CDF 
data, incorporating both direct and indirect production through the CSM (as a "fixed" 
contribution) which, in fact, is dominant at low and even moderate pt- 



3.1 Extraction of the colour-octet MEs 

In order to assess the effect of AP evolution on the fit parameters we show in table 1 
two sets of numerical values for the relevant colour-octet MEs obtained from a best x 2 
fit to Tevatron data |I2j using the CTEQ2L PDF: (i) the first row corresponds to a 



generation without AP evolution; (ii) the second set does take it into account. Notice the 
increase of < 8 T(1S) ( 3 5i) > in the latter case w.r.t. AP off (but to a lesser extent than 
for charmonium [|14j]) whereas M^ 15 ^ decreases consequently to keep the fit at low and 
moderate px values. 

Let us stress that our MEs numerical estimates have to be viewed with some cau- 
tion because of the theoretical and "technical" (due to the Monte Carlo assumptions) 
uncertainties. For example our algorithm for AP evolution should be regarded as a way 
to reasonably steepening the high-p T tail of the (leading-order) differential cross section 
which otherwise would fall off too slowly as a function of pt- 
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Table 1: Colour-octet matrix elements (in units of 10 -3 GeV 3 ) from the best fit to CDF data 
at the Tevatron on prompt T(1S) production. The CTEQ2L PDF was used with initial-state 
radiation on, and AP evolution off and on respectively. For comparison we quote the values 



given in [16, 11 1: 480 and 40±60 respectively. 



ME: 


< 8 T(15) ( 3 Si) > 


^T(lS) _ 5% V<O s I(li,, ( 3 P )> <O 8 I(li ' , ( 1 S )>\ 
5 y ml 5 J 


AP off 


93±18 


17±20 


AP on 


139±31 


6±18 



In figure 1 we show the theoretical curves obtained from our fit to CDF data (inde- 
pendently with AP evolution off and on in the generation) for both colour-singlet and 
colour-octet contributions. Let us remark that due to the px cut-off parameter set in the 
generation, only those experimental points with px > 1 GeV were used in the fit. Very 
good fits, with \ 2 /Ndf values close to unity, were found. 



3.1.1 Separated production sources for px > 8 GeV 

Current statistics does not permit to subtract indirect production sources to obtain the 
direct T(IS') production cross section along the full accessible p^-range. Nevertheless, 
feeddown from higher states (T(nS), Xbj{nP)) was experimentally separated out for px > 
8 GeV [|H| . We use this information to check our analysis a posteriori (rather than using 
it as a constraint in the generation) and to draw some important conclusions. To this end 
the relative fractions of the contributing channels for px > 8 GeV are reproduced in table 
2 from Ref. ||13|| . On the other hand we show in table 3 the fractions found in this work 



corresponding to the different generated channels for px > 8 GeV, following the notation 
introduced in section 2.1. 



Table 2: Relative fractions (in %) of the different contributions to T(1S) production from CDF 



data at px > 8 GeV [13]. Statistical and systematic errors have been summed quadratically. 



contribution 


Tevatron results 


direct T(1S) 


51.8±11.4 


T(2S)+T(3S) 


10.7±6.4 


Xb(lP) 


26.7±8.1 


Xb(2P) 


10.8±4.6 



By comparison between tables 2 and 3 we can conclude that the T(1S) indirect pro- 
duction from Xb/s decays is almost completely accounted for by the CSM according to the 
assumptions and values of the parameters presented in Section 2. Indeed, experimentally 
37.5±9.3% of T(IS') production is due to x w (lP) and x 6 j(2P) decays [R| while from our 
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Table 3: Relative fractions (in %) of the different contributions to T(1S) production at the 
Tevatron for px > 8 GeV from our generation. Possible contributions from Xbj{3P) states were 
not generated. 



contribution 



our generation 



T(1S) 



35^ 



42.3 



T(1S) I 



3(8) 



,(8) 



3.7 



T(1S) 



CSM 



14.9 



T(2g)+T(3g)| 

Xb(lP)\ CS M 

Xb(2P)\ CSM 



CSM 



3.0 
21.4 



generation we find a close value, namely 36.1%, coming exclusively from colour-singlet 
production! Moreover, assuming that a 7.7% from the 42.3% fraction corresponding to 
the colour-octet 3 S[ 8 ^ contribution (as expressed in Eq. (6)) can be attributed to the 
T(2S) + T(3S) channel in addition to the colour-singlet contribution (3%), we obviously 
get the fraction 10.7% for the latter, bringing our theoretical result into agreement with 
the experimental value. Furthermore, this single assignment implies to reproduce very well 
the experimental fraction (« 52%) of direct Y(IS') production by adding the remaining 

3 S[ 8 ^ contribution to the T(lS , )| ls ,(8) +3p (8) and T(1S) CSM channels (~ 53%). 

Of course all the above counting was based on mean values from table 2 and subject 
to uncertainties. Nevertheless, apart from the consistency of our generation w.r.t. experi- 
mental results under minimal assumptions, we can conclude, as an important consequence, 
that there is almost no need for Y(IS') indirect production from feeddown of Xbj states 
produced through the colour- octet mechanism. In other words, the relative contribution 
from P-wave states to < Oj^ 1S \ 3 Si) > \ tot in Eq. (6) should be quite smaller than naively 
expected from NRQCD scaling rules compared to the charmonium sector, in agreement 



with some remarks made in ||16|| . The underlying reason for this discrepancy w.r.t other 
analyses fill] can be traced back to the dominant colour-singlet contribution to the cross 
section at pr values as much large as 14 GeV (see figure 1) caused by the effective kr 
smearing - already applied to charmonium hadroproduction by one of us @] . 

On the other hand the corresponding velocity scaling rule in the bottomonium sector 
is nicely verified. Indeed from the value < ( 3 Si) > \ tot = 11.1 GeV 3 and the result 

found for < Oj^ 3 ^) > L, = 0.139 GeV 3 shown in table 1, the ratio 



<Ol {1S \ 3 S 1 )> 



tot 



0.012 (11) 



tot 



is in accordance with the expected order of magnitude ~ v 4 , where v is the relative velocity 
of the bottom quark inside bottomonium (v 2 ~ 0.1). 
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Figure 2: left : Predicted prompt T(1S) differential cross section at the LHC using the CTEQ2L 
PDF and AP evolution incorporated in the generation. A rapidity cut \y\ < 2.5 was required 
for bottomonium; dot-dashed line: contribution. Solid line: all contributions, right : 

Integrated cross section. 

4 T(15) Production at the LHC 

Bottomonium hadroproduction is especially interesting to check the validity of the colour- 
octet model as often emphasized in the literature HI. This becomes particularly clear 



at the LHC since experimental data will spread over a wider p^-range than at the Tevatron. 

Keeping this interest in mind, we generated prompt T(IS') resonances in proton-proton 
collisions at a center-of-mass energy of 14 TeV by means of our code implemented in 
PYTHIA employing the same colour-octet MEs of table 1 with AP evolution on. We 
present in figure 2 our theoretical curves for the T(IS') differential and integrated cross 
sections as a function of pt, including both direct production and feeddown from higher 
resonance states. 

In figure 3 we show our prediction for direct Y(IS') production. This is especially 
interesting if LHC detectors will be able to discriminate among such different sources of 
resonance production. To this end we generated T(IS') events through both the CSM and 
COM making use of the following parameters 

• <0^ 1S \ 3 S 1 ) 



> 



direct 



< o 



T(lS)/3 



Si) > 



direct 



9.28 GeV 3 (from | 
0.114 GeV 3 



M, 



T(15) 



6.0 GeV 3 



The first value corresponds to the CSM ME for direct production. The < Oj^ 1S \ 3 S\) > 
ME was obtained after removing the T(2S) + T(35) contribution according to the dis- 
cussion made in section 3.1.1, i.e. under the assumption that a fraction 7.7% from the 
42.3% in table 3 should be assigned to indirect production. Finally the M^ 1S ^ value is 
based on the assumption that this channel mainly contributes to direct production. 
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Figure 3: The same as in figure 2 for direct T(15) production at the LHC. 



5 Summary 

In this paper we have analyzed CDF measurements on T(IS') production cross sections 
at the Tevatron in a Monte Carlo framework. Higher-order QCD effects such as initial- 
state radiation of gluons and AP evolution of splitting gluons into (bb) states were taken 
into account. On the other hand, since different sources of T(1S) production were not 
experimentally separated along the full accessible p^-range we have included all of them 
in the generation and later fit. Only for p^ > 8 GeV feeddown from Xbj states was ex- 
perimentally separated out from direct production. We used such results as a consistency 
check of our analysis and to draw some conclusions summarized below. 

The numerical value of the < Oj^ 1S \ 3 S\) > \ tot matrix element should be ascribed 
almost totally to T(nS) states. This finding may be surprising when confronted with 
other results obtained from previous analyses pT| , IB| , where the contribution to the 



T(IS') yield through the colour-octet Xbj channels was thought as dominant [jT^, |TT| . 
On the contrary, we concluded from tables 2 and 3 that the colour- singlet production 
by itself can account for the feeddown of T(IS') from Xbj states. (Notice however that 
experimental uncertainties still leave some room for a possible COM contribution but to 



a much lesser extent than previously foreseen [H, |T5|].) On the other hand the different 
production channels are consistent (or can be made consistent) with the experimental 
relative fractions shown in table 2, after some reasonable assumptions. 

We have extended our study to LHC collider experiments (y^i = 14 TeV center-of- 
mass energy). In figure 2 we present our predictions for prompt production rates (i.e. 
including direct and indirect production) while in figure 3 we show our prediction for 
direct production alone. 

Lastly we conclude that the foreseen yield of T(lS')'s at LHC energy will be large 
enough even at high-p r to perform a detailed analysis of the colour-octet production 
mechanism and should be included in the B-physics programme of the LHC experiments, 
probably deserving (together with charmonia) a dedicated data-taking trigger. 
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